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From soft-collinear effective theory one can derive a factorization formula for the e+e^ thrust 
distribution da/dT with T = 1 — T that is appHcable for all T. The formula accommodates avail- 
able ff{(X^) fixed-order QCD results, resummation of logarithms at N^LL order, a universal non- 
perturbative soft function for hadronization effects, factorization of nonperturbative effects in 
subleading power contributions, bottom mass effects and QED corrections. We emphasize that 
the use of Monte Carlos to estimate hadronization effects is not compatible with high-precision, 
high-order analyses. We present a global analysis of all available e+e^ thrust data measured at 
2 = 35 to 207 GeV in the tail region, where a two-parameter fit can be carried out for as{mz) 
and ill, the first moment of the soft function. To obtain small theoretical errors it is essential 
to define i2i in a short-distance scheme, free of an ^(Aqcd) renormalon ambiguity. We find 
a,(mz) = 0.1135 ± (0.0002)expt± (0.0005)ai ± (0.0009)pert with x^/dof = 0.9. 
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A traditional method for testing the theory of strong interactions (QCD) and to make precise 
determinations of the strong coupling is the analysis of event-shapes measured at e+ collid- 
ers One of the most frequently studied event-shape variables is thrust 

T = max{ , (1) 

LilPil 

where the sum / is over all final-state hadrons with momenta pi, and the unit vector t that maximizes 
the RHS of Eq. ([T]) defines the thrust axis. It is convenient to use the variable T = I — T. For the 
production of a pair of massless quarks at tree level da/dx <x 5{x), so the measured distribution 
for T > involves gluon radiation and is highly sensitive to the value of a^- For t values close to 
zero the event has two narrow pencil-like, back-to-back jets, carrying about half the center-of-mass 
(cm.) energy into each of the two hemispheres defined by the plane orthogonal to t. For T close 
to the kinematic endpoint 0.5, the event has an isotropic multi-particle final state containing a large 
number of low-energy jets. The thrust distribution can be divided into three regions, 

peak region: t ~ 2Aqcd/ Q , 
tail region: lAQco/e < T < 1/3 , 
far-tail region: 1/3<T< 1/2. 

For T < 1/3 the dynamics is governed by three different scales. The hard scale Hh — Q, set by the 
e^e^ cm. energy Q, the jet scale, IJLj ^ Q\pi, the typical momentum transverse to t of the particles 
within each of the two hemispheres, and the soft scale \ls — Ql, the typical energy of soft radiation 
between the hard jets. In the peak region the distribution shows a strongly peaked maximum. 
Since T ^ 1 one needs to sum large (double) logarithms, (a/ln*T)/T, and do jdx is affected at 
leading order by a nonperturbative distribution, called soft function 5™°''. In the analysis presented 
in this talk we consider the tail region. It is populated predominantly by broader dijets and 3-jet 
events. Here the three scales are still well separated and one still needs to sum logarithms, but now 
/I5 S> Aqcd so soft radiation can be described by perturbation theory and the first moment of the soft 
function Q.\ = / d^(^/2)5'™°'*(/: — 2A). Many previous event-shape analyses have relied on Monte- 
Carlo (MC) generators to quantify the size of nonperturbative corrections. This is problematic since 
the partonic contributions implemented in MC generators are (i) based on LL parton showers and 
(ii) contain an infrared cut below which the perturbative parton shower is switched off and replaced 
by hadronization models that are not derived from QCD. Thus MC hadronization effects are not 
compatible with high-order perturbative event-shape predictions using the common MS scheme. 

In this talk we present a new analysis of e^e~ thrust data using the soft-coUinear effective 
theory (SCET), an effective theory for jets [|3|], to derive the theoretical QCD prediction of the 
thrust distribution. Within SCET it is possible to formulate a factorization theorem that allows to 
describe the thrust distribution for all T. The formula we use is iQ] : 

Due to lack of space we describe in the following only the main features of Eq. (^). For details, 
explicit analytic expressions, how our implementation improves upon earlier analyses in the liter- 
ature, and a complete set of references we refer the reader to Ref. iQ]. The term dds/dr contains 



2 



Global Fit ofasijnz) to Thrust at NNNLL Order with Power Corrections 



Andre Hoang 





cusp 


non-cusp 


matching /?[cts 


nonsingular 


7 a'" 


5 


LL 


1 




tree 


1 








NLL 


2 


1 


tree 


2 








NNLL 


3 


2 


1 


3 


1 


1 


1 


N^LL 




3 


2 


4 


2 


2 


2 


NLL' 


2 


1 


1 


2 


1 


1 


1 


NNLL' 


3 


2 


2 


3 


2 


2 


2 


N^LL' 


4Padc 


3 


3 


4 


3 


3 


3 



order 



as{mz) (No Gap) a5(mz) (With Gap) 



NLL' 


0.1203 ±0.0079 


0.1191 ±0.0089 


NNLL 


0.1222 ±0.0097 


0.1192 ±0.0060 


NNLL' 


0.1161 ±0.0038 


0.1143 ±0.0022 


N^LL 


0.1165 ±0.0046 


0.1143 ±0.0022 


N'^LL' 


0.1146 ± 0.0021 


0.1135 ±0.0009 


N='LL'(no qed) 


0.1153 ±0.0022 


0.1141 ±0.0009 


N^LL'(no o-fc/qed) 


0.1152 ±0.0021 


0.1140 ±0.0008 



Figure 1: (a) Ingredients for primed and unprimed orders used in our analysis. The numbers give the 
loop orders for the cusp and non-cusp anomalous dimensions, matching/matrix element contributions, the 
a^-running, the nonsingular distribution, the gap-anomalous dimensions, and the perturbative 7?-scheme 
subtractions 5 for our scheme for The 4-loop cusp anomalous dimension required at N^LL' order is 
estimated from Pade approximants. The associated uncertainty is negligible, (b) Central values and theory 
uncertainties for the fits at the different orders with and without the gap and renormalon subtractions. 



the singular partonic contributions. It factorizes further into a hard coefficient, a jet function 
and a partonic soft function governed by the renormalization scales Hh, pLj and /i^, respectively, 
and renormalization group (RG) evolution factors that sum logarithms between the hard, jet and 
soft scales. Using results from the existing literature, SCET allows to sum the logarithms at N^LL 
order [Q], which is two orders beyond the classic resummation method that is valid up to NLL 
order. The jet and partonic soft functions contain a/ [ln'^(T)/T]+ and a, 5(t) distribution terms. 
They are known to ff{aj), and at ^{cc^) all logarithmic terms are known from the renormalization 
group. Two unknown &{(X^) non-logarithmic constants contribute to the theory error in our high- 
est order numerical analysis. The hard function in our analysis is fully known at &{(X^) [14, 15] 
and also includes the axial- vector singlet contributions at ff{aj). To achieve a definition of the 
soft function moment Hi that is free of a Aqcd renormalon ambiguity, dds/dr contains subtrac- 
tions that eliminate partonic low-momentum contributions [^, ^. This requires the introduction 
of the additional scale-dependent model parameter A(/Xs) (with /i^ ~ jXs), called the gap param- 
eter, visible in Eq. @. In our numerical tail-data fits A(/Xs) is contained in Q.i. The evolution of 
A.{ij.r) follows a new type of infrared RG equation formulated in Refs. [^]. We have also included 
final-state QED matrix elements and QED RG corrections at NNLL order, derived from the QCD 
results. The term ddns/dT, called the nonsingular partonic distribution, contains the thrust dis- 
tribution in strict fixed-order expansion up to &{(X^) with the singular terms contained in dds/dr 
subtracted to avoid double counting. At ff{as) the nonsingular distribution is known analytically, 
and at ^{ocf) and ^{a^) we rely on numerical results obtained from the programs EVENT2 [ [l0| ] 
and EERAD3 [jll]] (see also [|l2|]). To achieve a consistent behavior in the far-tail region infrared 
subtractions need to be implemented here as well. A list of the perturbative ingredients for the dif- 
ferent orders we consider is given in Tab. |I|a. N^LL' is the highest order we consider and contains 
all currently available perturbative information. Finally, Ada^ / dr contains corrections to the singu- 
lar and nonsingular distributions due to the finite Z7-quark mass, using Refs. [13] for the consistent 
treatment and resummation for the singular terms. The entire partonic distribution is convoluted 
with the soft function S^°^ that describes the nonperturbative effects coming from large-angle soft 
radiation and can be determined from experimental data. The last term in the brackets indicates the 
parametric size of the dominant power corrections not contained in the factorization formula. For a 
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Figure 2: Plots of Hi vs as(mz)- (a) Includes perturbation theory, resummation of the logs, the soft model 
function and £2i with renormalon subtractions at jxr = 2 GeV. (b) As (a) but in a scheme Q.i without a gap, 
which gives perturbative results without the corresponding renormalon subtractions. The shaded regions 
indicate the theory errors at NLL' (brown), NNLL (magenta), NNLL' (green), N^LL (blue), N^LL' (red). 
The dark red elUpses in (a) and (b) represent the (x^i^ error ellipses for the combined theoretical, 
experimental and hadronization uncertainties. The ellipse in (a) is displayed again in Fig. 3b. The best fit 
points at N^'LL' with gap and renormalon subtractions shown in red in (a) each have X^/dof ~ 0.90. 



proper summation of large logarithmic terms it is necessary to adopt T-dependent profile functions 
for the renormalizations scales Hh, l^j, pi-s and jlR that follow the scaling arguments given above. 
For T — ;> 0.5 all profile functions need to merge into the hard scale }Xh to ensure that in the large-T 
endpoint region the partonic distribution coincides with the fixed-order result, so that it does not vi- 
olate the proper behavior at multi-jet thresholds. The variations of these profile functions estimate 
higher order perturbative uncertainty, and constitute our major source of theory uncertainty. 

In our analysis we fit the factorization formula (^ in the tail region to all available e^e^ 
thrust data from cm. energies Q between 35 and 207 GeV. In the tail region the distribution can 
be expanded in Aqcd/(2t) and thus described to high precision using as{mz) and Q.^. We carry 
out a two-parameter fit for these two variables. Fitting for D.\ accounts for hadronization effects in 
a model-independent way. For the fitting procedure we use a ;f ^-analysis, where we combine the 
statistical and the systematical experimental errors into the correlation matrix, treating the statistical 
errors as independent. We also account for experimental correlations of thrust bins obtained at one 
Q value by one experiment through the minimal overlap model, and find a similar central value 
to a completely uncorrelated treatment. To estimate the theoretical errors in the as — O-y plane 
we carry out independent fits for 500 different sets of theory parameters (for two unknown &{(X^) 
non-logarithmic constants, the four-loop cusp anomalous dimension, numerical uncertainties for 

2 3 

the &{ar) nonsingular distributions, parameters of the profile functions/renormalization scales) 
which are randomly chosen in their natural ranges with a flat distribution. We take the area covered 
by the points of the best fits in the a, — ^i plane as the theory uncertainty. 

The result of our fits for our default thrust tail range 6/2 < T < 0.33 (487 bins), at the five 
different orders we consider is displayed in Fig. ^. The left panel shows the results including the 
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Figure 3: (a) iXi-aim + l)-ellipse for the central fit at N^LL' order obtained from the experimental correlation 
matrix with default values for the theory scan parameters, (b) {Xmm + l)-ellipses of central fits for many 
different T-ranges in the tail region. The exponents display the number of data bins for each fit. The big red 
ellipse is the combined experimental and theoretical error ellipse that should be understood as 1-sigma for 
a,.(mz). 

gap and renormalon subtractions and the right panel without the gap and renormalon subtractions. 
Each dot corresponds to a best fit for a given set of theory parameters. The shaded areas envelop 
the best fit points and give the theory uncertainties. The numbers for central values and theory 
errors at each order are collected in Tab. [I]b and also display the size of the QED and Z^-quark mass 
effects. We see the excellent convergence of the fit results and the decrease of the respective theory 
uncertainties with increasing perturbative order. Moreover, including the gap and the renormalon 
subtractions leads to uncertainties that are about a factor of two smaller at the highest three orders. 
This illustrates the impact of the renormalon contributions and the necessity to subtract them from 
the partonic distribution. Our scan method is more conservative than the traditional error-band 
method. It is also important to quantify the experimental uncertainties. In Fig. ||a the {X^^ + 1)- 
ellipse for the central best fit at N^LL' order is displayed. For as{mz) we get a purely experimental 
error of (5as.)exp = 0.0002 and an error from the variations of D.i of {5as)cl^ = 0.0005. The 
latter uncertainty represents the hadronization error. Thus the theoretical uncertainties are about 
twice the hadronization error and about 4 times larger than the combined statistical and (correlated) 
systematic experimental errors. The dark red "circle" shown in Fig. ^ represents the total error 
including experimental, theoretical and hadronization errors. We note that to obtain stable fit results 
in the a, — Q.i plane it is essential to simultaneously fit data from different cm. energies Q because 
there is a strong theoretical degeneracy between 01^ and Hi. For each Q value an increase of can 
be compensated for the thrust distribution by a decrease of Hi. The strength of the degeneracy has, 
however, a strong dependence on Q, and can therefore be lifted by considering data from many 
different Q values within a single global fit. 

Finally, let us have a look at the dependence of the fits on the T-ranges used for experimental 
data. In Fig. |3|b the central best fits and the corresponding {x^in + 1) error ellipses for various 
T-ranges are displayed. The distribution of central fits and the ellipses is a remnant of the as — D-i 
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degeneracy just mentioned and arises from the dependence of how the degeneracy is lifted in a 
global fit on the selected T fit range. Including more peak data at small T leads to smaller 
(but sensitivity to the second moment Q.2 grows), and including less data increases the experimen- 
tal/hadronization errors. The distribution of the different best fit points represents a theoretical 
uncertainty which should not be double-counted with the theory uncertainty we already estimated 
from the parameter scan shown in Fig. This is compatible with the combined theoretical, ex- 
perimental and hadronization error from our default T-range also shown in Fig. |3|b. Our final result 
from our global analysis reads 

a,{mz) = 0.1135 ± (0.0002)expt ± (0.0005)^, ± (0.0009)pert 

= 0.1135 ± (0.001 l)tot. (3) 
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